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The Identification of EGRET Sources with Flat-Spectrum Radio Sources 
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ABSTRACT 



We present a method to assess the reliability of the identification of EGRET sources with 
extragalactic radio sources. We verify that EGRET is detecting the blazar class of AGN. 
However, many published identifications are found to be questionable. We provide a table of 
42 blazars which we expect to be robust identifications of EGRET sources. This includes one 
previously unidentified EGRET source, the lensed AGN PKS 1830—210 near the direction of the 
Galactic center. We provide the best available positions for 16 more radio sources which are also 
potential identifications for previously unidentified EGRET sources. All high Galactic latitude 
EGRET sources > 3°) which demonstrate significant variability can be identified with flat 
spectrum radio sources. This suggests that EGRET is not detecting any type of AGN other 
than blazars. This identification method has been used to establish with 99.998% confidence 
that the peak 7-ray flux of a blazar is correlated with its average 5 GHz radio flux. An even 
better correlation is seen between 7-ray flux and the 2.29 GHz flux density measured with VLBI 
at the base of the radio jet. Also, using high confidence identifications, we find that the radio 
sources identified with EGRET sources have larger correlated VLBI flux densities than the 
parent population of flat radio spectrum sources. 



Subject headings: Gamma Rays: observations — galaxies: active — quasars: general 



1. Introduction 

The Energetic Gamma Ray Experiment Telescope (EGRET) aboard the Compton Gamma Ray 
Observatory is sensitive in the energy range 30 MeV to 30 GeV (Thompson et al. 1993). A catalog of 
EGRET sources based on the first 30 months of exposure is given by Thompson et al. (1995, henceforth 
TH95). It includes 129 sources. Five are pulsars, 40 are listed as "high-confidence identifications" of AGN, 
and 11 are listed as "lower-confidence identifications" of AGN. A total of 71 are unidentified. 



^Astronomy Department, Boston University, 725 Commonwealth Ave., Boston, MA 02215 

^Former affiliations: Astronomy Department, University Of Maryland; Universities Space Research Association; Compton 
Observatory Science Support Center, NASA/GSFC; Gamma-Ray Astrophysics Branch, NASA/GSFC 

^Center for Astrophysics, Cambridge, MA 02138 

"'Dept. Of Phys. & Astron, University of Iowa, Iowa City, lA 52242 

SCode 662, NASA/GSFC, Greenbelt MD 20771 

^Max-Planck Institut fiir Radioastronomie, Bonn, GERMANY 



- 2 - 



We demonstrate quantitatively in this paper that the claim that EGRET is detecting the 7-ray emission 
of some members of the blazar class of AGN is valid (see von Montigny et al. 1995 for a compendium of 

EGRET blazar results for the first ^1/2 of the EGRET mission). By blazars, wc mean AGN with strong, 
compact, flat-spectrum (a be ^ —0.5, where S^u) oc j/") radio emission; and which also show continuum 
domination of the optical emission, and/or significant optical polarization, and/or significant changes in 
optical flux on short time scales (designated Optically Violently Variable, or OVV). These observational 
properties are believed to result from emission by material in a relativistic jet which is directed within 
^^10° of the line of sight. The blazar class includes objects classified as BL Lacertae type objects, high 
polarization quasars (HPQ), and OW quasars. The apparent 7-ray luminosity is as much as one hundred 
times larger than that at all other wavelengths for some flaring EGRET blazars. Variability of the 7-ray 
flux from some blazars on a time-scale as short as 4 hours (Mattox et al. 1997) implies that the 7-rays are 
emitted from a compact region. The 7-ray emission is not understood. It may offer fresh insight into the 
blazar phenomenon, and useful diagnostics. 

A careful analysis of the EGRET source identifications is in order because the source of the 7-ray 

emission is not well located. EGRET position estimates are imprecise because of the wide point spread 
function of the EGRET instrument. The half-angle of a cone which contains 68% of the EGRET events 
from a point source at a specific energy is well fit (Mattox et al. 1996a) by 

6»68 = 5.85°[£^/100MeV]-°-^3'' (1) 

The localization expected for a source is approximately 9 ~ 9psf/VN where N is the number of 7-rays 
detected from the source (Thompson 1986). For weak EGRET sources, the 95% confidence position error 
ellipses are as large as ^5 square degrees. 

In the first EGRET catalog (Fichtel et al. 1994), EGRET blazars were classified as "positive detections" 
or "marginal detections" depending on whether the significance of point-source 7-ray emission exceeded 5a 
or 4(7 respectively. The confidence of the identification was not explicitly addressed. The second EGRET 
catalog (TH95) did explicitly address the confidence of identification. TH95 classified an identification as 
"high-confidence" if the radio source was located within the 95% confidence contour of the EGRET position 
estimate. Also, TH95 classified some radio sources beyond the 95% confidence contour as "lower-confidence" 
identifications. The 5 GHz flux density (5*5) of a radio counterpart was generally required by TH95 to be 
^ 1 Jy; and the radio spectral index to be generally a ^ —0.5. No formal consideration was given to the 
size of the EGRET error region and where the radio source lay within that error region (other than whether 
the source was located within the 95% confidence contour of the EGRET position estimate region or not). 
No formal consideration was given to the number density of potentially confusing radio sources (other than 
a loose requirement that iSs be ^1 Jy). No formal consideration was given to the radio spectral index 
(other than a loose requirement that a be ^ —0.5). No consideration was given to the a priori probability 
of detecting a radio source by EGRET it is less than unity because many blazars which are bright at 
other wavelengths have been observed by EGRET without being detected. Also, the decision by TH95 on 
which radio sources beyond the 95% confidence contour to include as "lower-confidence" identifications was 
subjective. 

We have developed an analysis of EGRET radio source identification which quantitatively incorporates 
all of this information. We avoid a cutoff in the search for an identification at an arbitrary confidence 
contour of the EGRET position estimate region. Also, we avoid terminating the search for candidates at 
a speciflc minimum radio flux density or radio spectral index. Bayes' theorem is used to incorporate the 
a priori probability of detecting a radio source (see Sturrock 1973 and Loredo 1990 for discussions of the 
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use of Bayesian statistics in astronomy). It is reasonable to use a Bayesian analysis here for the following 
reasons: 1) we can demonstrate with very high confidence (see §3.1) that EGRET is detecting flat-spectrum 

radio sources without making any a priori assumptions; 2) wc can derive a reasonable estimate for the a 
priori probability of EGRET detecting a radio source as a fmiction of radio flux (sec §3.2). 

The method is described in §2. In §3 we demonstrate conclusively that EGRET is detecting blazars 
and obtain an estimate for the a priori probability of EGRET detecting a radio source as a function of radio 
flux; and then provide a quantitative study of the reliability of the EGRET identifications of TII95 and 
Thompson et al. (1996), and search for additional radio identifications among the previously imidentified 
EGRET sources. Limitations of this analysis are discussed in §3.6. Implications are discussed in §4. In 
§3.3.1, we use the high probability identifications of §3.3 to demonstrate that the EGRET blazars have 
higher VLBI flux densities than the parent population of flat-spectrum radio sources. 

2. EGRET Source Identification 

We analyze the probability of a correct identification of an EGRET source with a radio source. The 
calculation considers the size of the EGRET error region and where the radio source lies within that error 
region. Baycs' theorem is used to include consideration of the number density of potentially confusing radio 
sources and the a priori probability of detecting a specific radio source by EGRET. This method is an 
extension of that used by de Ruiter, Willis and Arp (1977) to identify optical counterparts of radio sources. 
This method has also been applied to identification of IRAS sources (Wolstencroft et al. 1992) and X-ray 
sources (Schachter, et al. 1997). 

2.1 Radio Catalogs 

Because blazars are flat-spectrum radio sources, high frequency surveys are appropriate for this study 
since they efficiently detect fiat spectrum sources. We have used the 4.85 GHz Greenbank (GB) survey 
(Condon, Broderick, & Seielstad 1989) for the northern sky and the 4.85 GHz Parkes-MIT-NRAO (PMN) 
survey (Griffith & Wright 1993) for the southern sky. Both surveys have a threshold flux density of ^30 
mJy and give positions with -^20" uncertainty. They are both confusion limited at Galactic latitudes less 
than 3°. Therefore, we do not attempt to identify the 19 unidentified EGRET sources at latitudes less 
than 3° in this paper. Also, the PMN survey is confusion limited near the Large Magellanic Cloud and 
Centaurus A. For this reason, EGRET sources 2EG J0532-6914 and 2EG J1324-4317 are also excluded from 
this analysis. 

The GB survey covered the declination range 0°< 6 <75°. We have used the catalog of Becker, White, 
and Edwards (1991) of 53522 sources derived from this survey. The White and Becker (1992) catalog of 
sources from the 1.4 GHz GB survey was used to obtain an a1.4_4.85 spectral index for the 22384 sources 
which are found in both catalogs. Because the 1.4 GHz survey preceded the 4.85 GHz survey by ~3 years, 
variability of blazar type AGN can result in this index differing by as much as -^1 from what would be 
obtained with simultaneous measurements. 

The PMN survey (Griflith & Wright 1993) covered the decHnation range -87.5°< S <10°. There are 
small gaps in the survey due to observational problems. The PMN catalog contains 49329 sources with 
55 > 30 mJy. Spectral indices are provided for the -^10% of the PMN sources which were detected in the 
Parkes 2.7 GHz survey (Bolton, Wright, and Savage 1979). Again, variability can cause a substantial error 
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in this spectral index. Only one unidentified source (2EG J1332+8821) falls completely outside of these 
surveys. The error eUipses of two EGRET sources (J1409-0742 and J1513-0857) fall partially in gaps 

in the PMN survey, hut both have strong identifications from the PKS catalog. We retain the original 
eight digit position names of both catalogs. The GB source names are prefixed with "B" because the 
position names are based on B1950 celestial coordinates. The PMN source names are prefixed with "J" 
corresponding to J2000 celestial coordinates. We use the GB survey in the region where the surveys overlap, 
0°< 6 <10°. 



2.2 Probability of Correct Identification 

For a radio source found at an angle r from the center of an EGRET error ellipse, there are two 

possibilities: the source is the radio counterpart of the EGRET source; or, the source is a confusing source 
which is close to the EGRET source by accident. We denote the probabilities of these two hypotheses by 
p{id\r) and p{c\r) respectively. 

If the radio source is merely a confusing source, the Poisson distribution pertains. The probability that 
one or more radio sources are accidently within the angle r is 

p(r|c) = l-e-'''/'^", (2) 

where ro is the characteristic angle between confusing sources, 

ro = {np{S5,a))-i, (3) 

and p{S5,a) is the number density of radio sources with at least the flux density of the candidate radio 
source, and a spectrum at least as flat. The differential probability is 

dp{r\c) =2^6-''^ /''odr (4) 

For sources without spectral indices, we assume a spectral index of —0.5. This results in a conservative 
estimate of the signiflcance of the identification for radio sources which are indeed flat spectrum. A refined 
estimate of the probability of correct identification can be made after further radio observations yield a 
spectral index. We obtain p{Ss, a) from the GB survey for the 16300 square degree region 0°< S < 75°, and 
|6| > 10°. The integral of p{S5,a) over a is observed to have a Euclidean slope between 100 and 1000 mJy 
in agreement with the result of Kellerman and Wall (1987). We use this number density also for analysis of 
PMN counterparts. 

If the radio source is the counterpart of an EGRET source, it is appropriate to use the EGRET data 

to estimate the probability distribution, p{r\id), of the angle r. A likelihood analysis of the EGRET data 
(Mattox et al. 1996a) yields the 95% confidence position contour from the likelihood ratio test using Wilks' 
theorem. These contours are usually fit accurately with ellipses. The ellipse is centered on the centroid 
of the 95% confidence location region which is given by TH95 as the position estimate. This position is 
observed to differ by only a small amount from the position of maximum likely which is formally the correct 
position estimate. The ellipse is specified by the semi-major axis a, semi-minor axis b, and the position angle 
of the major axis ^. The departure from circularity is often substantial. This results from source confusion 
(nearby 7-ray sources below the EGRET detection threshold) . error in the model of Galactic diffuse 7-ray 
emission, and statistical fluctuation. The elliptical fits given by TH95 were used in our analysis. The 
parameters for the elliptical fits which are not given in TII95 (because of an imprecise fit) have also been 
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used in this analysis: 2EG J0426+6618, a =66', b =54', (j) =27°; 2EG J0511+5523, a =135', b =36', <j) =24°; 
2EG J0545+3943, a =48', b =47', </> =94°; These fits are all accurate to within --30% in radius which is not 
substantially worse than the potential effect of error in the Galactic diffuse model and undetected nearby 
7-ray sources. 

For strong sources, we observe that the dependence of the log of the likelihood of the EGRET data on 
the assumed source position is well represented by a paraboloid. The 95% confidence contour corresponds 
to a drop in the log of the likelihood of 3.00 from the maximum (Mattox et al. 1996a). Therefore, we 
represent the log likelihood surface corresponding to an elliptical error region as an elliptic paraboloid 
In Ljnax — Ini = 3r^/^'^, where 

= 0082(6* -(/))+ b-^ sin^id - (/))]-^/^ (5) 

and 9 is the position angle. For a circular error region, 5* is simply the 95% confidence radius. From 
Wilks' theorem, 2{lnLmax — InL) is distributed as in the null hypothesis (Mattox et al. 1996a). Thus, 
the confidence with which a source may be expected within a region delineated by a A = In L^ax — In L 
decrease in log likelihood is 

/>oo 

C=l- X2(C)dC = l-e-^ (6) 
Therefore, the probability of the EGRET source being located at an angle smaller than r is 

p{r\id) = 1 - e-^-V*' (7) 

and 

dp(r|id) =6^e^3'''/*'dr (8) 

Bayes' theorem allows the observed angle between the radio source and the EGRET position estimate, 
r, to be used with an a priori probability and the expected distributions of r under the two possible 
hypotheses (equations ^ and ^ to obtain p{id\r), the a posteriori probability that a radio source is the 
correct identification of an EGRET source. 

p{id\r) = ^ p{id)dp{r\id) 

dp{r) p{id)dp(r\id) + p(c)dp(r\c) 

With the substitution rj = p{id) for the a priori probability, and using the fact that one of the two 
hypotheses is true, either a source is the identification or it is a confusing source, p{c) = 1 —p{id); we obtain 

-^LR 

pHO = , . ' (10) 



1-1] 



where the likelihood ratio, LR, is defined to be 



dp{r\id) _ rg .(3^^._,^,.) 

dpir\c) " ^ ' 

Figure 1 shows the LR and rj dependence oi p{id\r) given by equation p^ . 

We note that LR diverges at large r if '5 > VS tq. In this case, the differential probability of r under 
the confusion hypothesis falls off more rapidly with r than it does under the identification hypothesis. This 
refiects the fact that the size of the EGRET error region, ^P, is larger than the characteristic angle between 
confusing sources, ro. A meaningful identification is not possible. We note that LR < 1 at r = in this 
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case. Therefore, we require < V3 vq before we consider a radio source as a candidate identification of 
an EGRET source. This Umits the number density of sources for which an attempted identification is 
meaningful. The consequence is that a 5 GHz flux density of ^ 50 mJy is required for a source to be 
a meaningful candidate for identification given the error ellipse sizes typical of the unidentified EGRET 
sources. 



3. Results 

Table 1 presents the results of the analysis described in §2 for the 102 EGRET sources from TH95 
which are not identified as pulsars, and which are not in regions where the radio surveys are source confused 

< 3°, and near the LMC and Gen A). All radio sources are given for which an a posteriori probability, 
p{id\r) > 0.001 is found. The radio source with the highest value of p{id\r) is included for all EGRET 
sources which fall within the radio survey regions. We note that sources with radio flux densities as small 
as ~50 mJy appear, but with very small probabilities of being EGRET sources. The same analysis is 
presented in Table 2 for the sources in the Supplement to the Second EGRET Catalog (Thompson et al. 
1996). The content of Tables 1& 2 is: 

Column (1) contains the J2000 position name of the EGRET source from TH95 or Thompson et al. 
(1996). 

Column (2) is the Galactic position of the EGRET source (from TH95 or Thompson et al. 1996). 
For the five sources where TH95 accidently gave the radio position (0528+134, Mkn 421, 3C 279, 
1406-076, 3C 454.3), the EGRET positions have been substituted. 

Column (3) is the variability index which is — Iog2o(-P), where P is the probability that the source 
is invariant based on a test of flux for all available viewing periods (McLaughlin et al. 1996). A 
larger variability index corresponds to a higher probability that the source is variable. 

Column (4) is the identification given (for some sources) by TH95. A question mark following the 
B1950 radio position name means that TH95 indicated that it was a "lower-confidence identification"; 
otherwise, the identification was indicated by TH95 to be a "high-confidence identification" . 

Column (5) is the name of the potential radio counterpart. The original B1950 position names are 
used for Greenbank sources, and the original J2000 position names are used for PMN sources. 

Column (6) is a more common name for some potential radio counterpart. 

Column (7) is the 4.85 GHz fiux density in mJy from the survey catalog. 

Column (8) is the radio spectral index, a, from the radio survey flux densities if available, S{v) oc v°' . 

Column (9) is ro, the characteristic angle (in arcmin) between confusing sources which are at least as 
bright and flat as this radio source (see equation ^) . 

Column (10) is 77 = p{id), the a priori probability that the radio source is a 7-ray source from equation 

H- 

Column (11) is the angle, r (in arcmin), between the EGRET position estimate and the radio position. 
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Column (12) is 5*, the radius (in arcmin) of the 95% confidence contour in the direction of the radio 
source from equation ^ 

Column (13) is the position confidence contour at the radio position (from equation^). 

Column (14) is the likelihood ratio (equation |ri|) indicating the strength of the indication for the 
identification. 

Column (15) is the a posteriori probability that the identification is correct as determined with 



3.1 The Statistical Significance of the EGRET Detection of Radio Loud AGN 

We can now address an important question: with what confidence can we conclude that EGRET is 
detecting radio loud AGN? At this point, it is not appropriate to use Bayesian statistics. We have not yet 
determined the a priori probability of a flat spectrum radio source being an EGRET source. As illustrated 
by Figure 1, the a posteriori probabilities depends very strongly on the a priori probability. Note that for a 
priori probabilities 77 = 1 or r/ = 0, equation |l^ shows that Bayes' theorem implies that the EGRET data 
do not influence the expectation about whether or not a radio source emits 7-rays. 

Since Bayesian statistics are of no use at this point, we apply a "frequentist" ansatz, a Monte Carlo 
calculation. We use LR (equation ^l]) as an empirical indication of the strength of a potential identification. 
We find that the distribution of LR for the most likely identification of each source in Table 1 differs 
dramatically from the distribution found by Monte Carlo calculation assuming no EGRET detections of 
radio loud AGN. Simulated EGRET source positions were chosen randomly with |6| > 3° and used with 
the error ellipse parameters of TH95. The analysis of §2 for ~3000 simulated sources yielded an estimate 
for the distribution of LR in the null hypothesis. A long tail at large LR was observed. Confusing radio 
sources with LR >1{)^ were found at a frequency of 0.010 ± 0.002. This corresponds to the occurrence at a 
non-negligible frequency of spurious coincidences which look fairly solid. We find no strong dependence of 
LR on the size of the error ellipse, supporting the validity of the method we describe in §2. 

The strongest identification (and coincidently the first, Hartman et al. 1992) is 3C 279 with 
LR — 4.6x10^. The large value of LR reflects the fact that 3C 279 is centrally located in a small error 
ellipse, and that sources as bright and fiat as 3C 279 are very rare: p{Sz > 17Jy,Q! > 0.3) = 4x10"'^ per 
square degree. The largest observed LR in the ~3000 Monte Carlo trials was 6x10'^. Therefore, the 3C 279 
detection alone indicates that AGN are seen by EGRET with a confidence well above (l-102/3000)=97%. 
With the admission of an a priori probability as small as rj =10"'^, equation |l^ implies an a posteriori 
probability of 0.998 that 3C 279 is the EGRET source counterpart; and with the self-consistent a priori 
probability of 77 = 0.2 which we derive in §3.2, the a posteriori probability that the identification is correct 



Even stronger evidence that EGRET is detecting radio loud AGN is obtained by considering all sources 
in Table 1. From the Monte Carlo calculation, LR > 100 occurs with a frequency of 0.05 ± 0.004 in the null 
hypothesis. In comparison, 32 of the 102 EGRET sources in Table 1 have an identification with LR > 100. 
The probability to obtain this by chance is 



equation hffl 



is 0.999991. 




(12) 
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This is a very conservative estimate of the significance because many of the values of LR are much larger 
than 100. However, there are "trials" made in the selection of the type of source to try to identify with 
EGRET sources. The number of such trials cannot exceed the ~100 independent types of astronomical 
objects. There are also possible "trials" made in the selection of the energy range (e.g. 2022—077) or diffuse 
model (e.g. 1622—253) which was used to obtain an EGRET position for an identification. This number 
is 'Z 2. The choice of the minimum value of LR of 100 for this analysis also corresponds to a number of 
"trials", but the number is less than ^ 10. Also, the use of the radio spectral index in the determination of 
the number density of radio sources can be regarded as a trials factor of size ~2. (We demonstrate in §4.2 
that EGRET is generally not detecting steep-spectrum radio sources, so the use of the radio spectral index 
in the determination of the number density of radio sources is justified.) The number of potential "trials", 
the product of these, is far less than the inverse of the result of equation |l^. Therefore, we conclude with 
very high confidence that EGRET is detecting radio loud AGN. We note that this result does not depend 
on any a priori assumption about the emission of 7-rays by AGN. 



3.2 The A Priori Probability of EGRET Detecting a Radio Source 

Now that we have clearly established that EGRET is detecting some radio-loud AGN, it is sensible 
to estimate the a priori probability, 77, for EGRET to detect a particular flat-spectrum radio source. This 



will allow the use of equation 10 to estimate the a posteriori probability of a correct identification of any 
specific source. We do not attempt here to model the effect on rj of the variation across the sky of the depth 
of the EGRET exposure. We expect that this cannot effectively be done because of the large uncertainty 
of the gamma-ray luminosity function of AGN (Chiang et al. 1995), but that it is not a huge effect. Also, 
we assume that 77 is not dependent on the radio spectral index, although we derive 77 for a > —0.5. It is 
appropriate to thus give sources with a < —0.5 the "benefit of the doubt" because of the potential error in 
a obtained from the surveys. 

We seek a functional dependence on radio flux for ry which is self-consistent as done by de Ruiter et al. 
(1977) in their identification of optical counterparts of radio sources. Specifically, we solve for a value of 
rj such that the integral of p{id\r) given by equation |l^ (for radio sources with a > —0.5) divided by the 
number of radio sources considered (with a > —0.5) yields the assumed value of rj. We analyze 77 for several 
S5 intervals. These intervals are chosen to provide adequate statistics while preserving sensitivity to the S5 
dependence of rj. 

For the flux density range ^5 > 5 Jy, iteration leads to the self-consistent value of r/ — 0.19 ± 0.08 
with 'Sp{id\r) = 6.0 for 32 radio sources. For the flux density range 2 < S'5 < 5 Jy, iteration leads to the 
self-consistent value of 77 = 0.13 ± 0.04 with T,p{id\r) — 14.7 for 109 radio sources. For the flux density range 
1 < S'5 < 2 Jy, iteration leads to the self-consistent value of 77 = 0.034 ± 0.011 with T,p{id\r) = 9.1 for 266 
radio sources. And, for the flux density range 0.6 < S'5 < 1 Jy, iteration leads to the self-consistent value 
of 77 = 0.0047 ± 0.0032 with T,p{id\r) = 2.2 for 473 radio sources. The uncertainty of 7/ is simply assumed 
to be r][Sp{id\r)]^^/^ . This is the approximate expectation for large S5 because of Poisson fluctuation. For 
small S5, it might be an overestimate, but we expect that it is reasonable given potential systematic error. 

We have chosen an analytical representation for 77 which has appropriate limiting values for both high 
and low values of S'5. The high limiting value is chosen to be consistent with 77 = 0.19 ± 0.08 for S5 > 
5 Jy. An additional two parameters have been fit to match the relatively well determined values of 77 for 
1 < S'5 < 2 Jy and 2 < S'5 < 5 Jy. 

77(S5) = 0.2(1 -e-0-"7(^^/Jy)"') (13) 
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This fit is shown along with the data in Figure 2. It is shghtly in excess of the data for 0.6 < 55 < 1 Jy, 
but only by a small amount relative to the uncertainty. It is not inappropriate to thus give small S5 sources 
the "benefit of the doubt" because of variability. We have also done this analysis for an alternate choice 
of 5*5 bins: 0.5 < < 0.8 Jy, 0.8 < < 1.5 Jy, 1.5 < S'5 < 4 Jy, & ^5 > 4 Jy. The values of 77 thus 
obtained are consistent with equation Also, the value of rj for S5 bins a factor of 2 narrower is consistent 
with equation Thus, the derived a priori probability is not strongly dependent on binning. In §4.5, we 
discuss the strong probability that the very low value of rj for 0.6 < 6*5 < 1 Jy implies a correlation between 
7-ray flux and radio flux. 

3.3 Verification of Previous Identifications 

With the method of §2 and the prescription for an a priori probability of equation we can now look 
objectively at the radio source identiflcations of TII95 which are noted in column 4 of Table 1. Most of the 
"high-confidence identifications" of TII95 are indeed so. Of these 40 identifications; six have p{id\r) > 0.999, 
seven have 0.99 < p{id\r) < 0.999, nine have 0.90 < p{id\r) < 0.99, ten have 0.50 < p{id\r) < 0.90, and six 
have 0.10 < p{id\r) < 0.50. Only two have p{id\r) < 0.10: 2356+196 and 1604+159. The 2EG JOOOO+2041 
EGRET 95% confidence position contour contains 3.2 square degrees and radio sources as flat and bright as 
PKS 2356+196 have a number density of 0.01 per square degree. We find that the fraction of flat-spectrum 
radio sources with the 5 GHz flux density of PKS 2356+196 (700 mJy) which are detected by EGRET is 
only about 0.006. Thus we flnd with 93% confidence that PKS 2356+196 is a confusing source rather than 
the identification of 2EG JOOOO+2041. It is certainly possible that 2EG JOOOO+2041 is PKS 2356+196, 
but we suggest that it should not be considered to be a "high-confidence identification" . Similarly, PKS 
1604+159 is also a weak 5 GHz sources in a large EGRET error ellipse. 

Four of the ten "lower-confidence identifications" of TH95 are found to be very unlikely with 
p{id\r) < 0.005. All are modest 5 GHz sources well beyond the 95% contour of large EGRET error ellipses. 
On the other hand, two of the ten "lower-confidence identifications" of TH95 are found to be quite likely. 
With p{id\r) > 0.95, the identifications of 0458-020 and 0234+285 are found to be better than 55% of the 
"high-confidence identifications". Both are bright, fiat-spectrum 5 GHz sources which have a low number 
density and the fraction of sources this bright at 5 GHz which are detected by EGRET is large. 

Although TH95 classify PKS 1406-076 as a "high-confidence identification", their method would also 
allow it to be classified as a "lower-confidence" identification because PKS 1406-076 lies right on the 95% 
confidence contour of TH95 (see the ApJ CD ROM). We find that the identification of PKS 1406-076 is 
correct with 86% confidence. Although PKS 1406-076 is not near the center of the error region, the 95% 
confidence position contour contains only 0.2 square degrees. The variability of the EGRET flux and the 
flnding of a possible correlation of the variability of optical flux with the 7-ray flux (Wagner et al. 1995) 
indicates strongly that the identiflcation is correct. 

Table 2 shows the analysis of §2 for the sources in the Supplement to the Second EGRET Catalog 
(Thompson et al. 1996). The three AGN identiflcations in the Supplement to the Second EGRET Catalog 
are found to be compelling. 

3.3.1 The VLBI Flux Density of EGRET Blazars 

We have examined the 2.29 GHz VLBI flux densities (Preston et al. 1985) for the sources in Tables 1&2 
which either have p{id\r) > 0.9, or have p{id\r) > 0.7 for a EGRET source which is highly variable (PKS 
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Fig. 1. — The relationship of equation (3) for the probabiHty of a correct identification p{id\r) as a function 
of the log of LR for the indicated values of 77. 
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Fig. 2. — The measured a 'priori probability r\ — p{id) of a radio source being detected by EGRET as a 
function of S5. The fit of equation [l^ is also shown. 
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1406—076) or for which the confidence of the identification is otherwise enhanced (Mrk 421). Thirty of these 
31 sources (all except PKS 1830-210) were measured by Preston et al. (1985). The 30 EGRET blazars 

are shown in comparison to the parent population in Figmc 3. The marginal distribution of the VLBI flux 
density is shown in Figure 4 for the sources in Figure 3 with a > —0.5. It is immediately apparent that 
the EGRET blazars have higher VLBI flux densities than the parent population. A Kolmogorov-Smirnov 
(KS) test shows that these distributions differ with 99.9994% confidence. VLBI upper limits were treated 
as detections in this KS test. This is conservative because Preston et al. (1985) obtained VLBI detections 
for all 30 of the EGRET blazars they observed. 

We note that the amount of VLBI flux density is a much better discriminant of EGRET blazars than 
the VLBI core visibility. Figure 5 shows the distribution of the visibilities from the catalog of Preston et 
al. (1985) for the EGRET blazars. A KS test indicates with only 81% confidence that the EGRET blazars 
have higher visibilities than the parent population. The VLBI visibility is influenced by the amount of flux 
from the radio lobes. If the properties of radio lobes do not reflect whether 7-rays are visible, it is expected 
that the variance of lobe radio luminosity will cause the observed effect. 

This is the first report of a radio property which distinguishes the EGRET blazars from the parent 
population of fiat-spectrum radio sources. It supports arguments that the 7-ray emission is taking place 
at the base of a relativistic jet. This finding also has a practical application. A large VLBI flux density 
enhances the confldence of an EGRET blazar identification. 



3.4 Identification Analysis of Previously Unidentified EGRET Sources 

An obvious hypothesis is that some of the unidentified EGRET sources are also blazar type AGN. 
Our results are shown in Table 1 for the unidentified EGRET sources in the Second EGRET Catalog 
(TH95) and in Table 2 for the unidentified EGRET sources in the Supplement to the Second EGRET 
Catalog (Thompson et al. 1996). We find evidence of additional blazar identifications. The strongest new 
identification is lensed blazar PKS 1830—210 near the direction of the Galactic center with p{id\r) = 0.98. 
A total of four new potential identifications with p{id\r) > 0.5 appear in Tables 1 & 2. They include the ~2 
Jy sources J1650— 5044 and J1802— 3940 which are within 10° of the Galactic plane and have not yet been 
sufficiently studied to be certain that they are not Galactic sources. These 4 likely identifications are given 
in Table 4 along with an additional 11 possible identifications with p{id\r) >0.05. Two additional sources 
are also included in Table 4 which have 0.02< p{id\r) <0.05 and an indication of variability with confidence 
in excess of 99% (McLaughlin et al. 1996) — variable EGRET sources are more likely to be blazars. 



3.5 Summary of Blazar Identifications 

Table 3 gives the 42 EGRET blazar identifications which we expect to be robust. All of these sources 
have been established to be AGN. Most of these identifications have p{id\r) > 0.9. Several sources are 
included with 0.3< p{id\r) <0.9 if there is high confidence of variability of the EGRET source (in excess of 
99%), or if a large VLBI flux density has been observed. Our minimum requirement is 300 mJy in the 2.29 
GHz survey of Preston et al. (1985) or in a subsequent VLBI observation (details given in the footnotes 
of Table 3). Prom Figure 4, the requirement of a 2.29 GHz VLBI flux density > 300 mJy excludes 54% 
of potential confusing sources. Table 3 includes 33 of the 40 "high-confldence identifications" of TH95, 
2 of 10 "lower-confidence identifications" of TH95, 1 new identification of an unidentified TH95 source, 
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Radio Spectral Index 



Fig. 3. — A scatter plot of the 2.29 GHz VLBI flux densities and radio spectral indices from Preston et al. 
(1985) for the 1259 of 1398 sources in their catalog with radio spectral indices. Upper limits for VLBI flux 
densities are shown for 430 of these. Those which arc robustly identified with EGRET sources {p{id\r) > 0.9, 
or p{id\r) > 0.7 for PKS 1406-076 and Mrk 421) are circled. The EGRET blazar shown with a = +1.3 is 
3C 454.3 — the index of Preston et al. (1985) probably differs from the GB survey value a = +0.1 because 
of variability. 




Correlated VLBI Flux Density at 2.29 GHz (mJy) 



Fig. 4. — The marginal distribution of the VLBI flux densities of Figure 3 for a > —0.5. The 882 sources 

in the catalog of Preston et al. (1985) with a > —0.5 are shown with the histogram using thin lines. Upper 
limits arc treated as detections (at the value of the upper limit). The peak at ~120 mJy results from this 
treatment of upper limits. The 30 EGRET blazars circled in Figure 3 are shown with the histogram using 
bold lines. 
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3 identifications from Thompson et al. 1996 (one of these, NRAO 190, was first identified by McGlynn 
et al. 1996), and 3 identifications from subsequent EGRET exposures. Table 4 comprises the potential 

identifications for which wc expect further investigation will be most useful. Four of these potential 
identifications are very compelling with p{id\r) > 0.5. Some sources warrant specific comments. 



B0115+0242: EGRET source 2EG J0119+0312 is identified with 6% confidence with radio source 
B0115+0242. Because there is a 92% confidence indication of variability for 2EG J0119+0312, this 
identification is actually more likely than 6%. The optical magnitude is 17.5. The red shift is 0.672. 
It shows OVV behavior. This source is also known as PKS 0115+02 and 3C 37. The a1.4_4.85 radio 
spectral index is -0.7. It would be worthwhile to conduct more radio and optical flux monitoring 
observations of this source to obtain a radio spectral index from simultaneous observations, to 
characterize variability, and to determine the amount of mm radio flux density. Also, it would be 
worthwhile to observe the parsec scale radio structure of this source, and determine a spectral index 
of possible compact radio emission. 

0219+4248 (3C 66A): We find p{id\r) of 4.0% for the identification of EGRET source J0220+4228 
with 3C 66A assuming a radio spectral index of —0.5. However, a radio index of ~0. is indicated 
by 22 and 37 GHz monitoring (Terasranta private communication). Using this radio spectral index, 
p{id\r) is found to increase from 4.0% to 6.2%. We note that 3C 66A is a BL Lac type blazar which 
shows optical variability and that the VLBI flux density at 2.29 GHz from Preston et al. (1985) is 
0.27+0.05 Jy. Thus, this identification is credible. However, Verbund et al. (1996) identify this source 
with binary ms pulsar J0218+4232 after finding ~2c7 evidence of modulation of the EGRET flux after 
epoch folding. 

0336-019: Source 0336-019 (CTA 26) was detected by EGRET during the interval 05/09/95 — 
05/23/95 (Mattox et al. 1995). The EGRET 95% confldence contour is well flt by an eUipse centered 
on 1=188.25, b= —42.34, and with semi-major axis o=44', semi-minor axis b= 37', and (f)= 27°. A 
complete description of multiwavelength observation is in preparation. 

0446+112: This identification by TH95 is plausible, but narrowly fails to meet our criteria with 
p{id\r) = 0.31 and Svlbi = 200 ± 20 mJy. It is possible that this 2.29 GHz VLBI flux density 
measurement (Preston et al. 1985) is under representative because of variability. 

1229—021: This identification by TH95 is plausible, but fails to meet our criteria with p{id\r) = 0.32 
and Svlbi = 90 + 30 mJy. It is possible that this 2.29 GHz VLBI flux density measurement (Preston 
et al. 1985) is under representative because of variability. 

B1324+2226: EGRET source 2EG J1324+2210 is identifled with 63% confidence with this GB 

source. The 0:1.4-4.85 radio spectral index indicated by the (non-simultaneous) GB survey flux 
densities is +1, very strongly inverted. Recent simultaneous cm radio measurements indicate a radio 
spectral index of ^0. (M. Aller, private communication). A VLBI position was published by Johnston 
et al. (1995) but a VLBI flux was not given. 

J1650-5044: Our analysis shows a strong probabflity (92%) that this PMN source is EGRET source 
2EG J1648— 5042. However, observations at other wavelengths are required to ascertain that this 
radio source is an AGN rather than a Galactic radio source. This source was previously detected as 
radio source MRC 1646—506. Otherwise, it does not appear in the hterature. An observation with 
the Australian Telescope Compact Array would be useful (it is too far south for the VLA). 
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J1802-3940: 2EG J1800-4005 is identified with PMN J1802-3940 with 81% confidence. 
Observations at other wavelengths are required to ascertain that this radio source is an AGN rather 
than a Galactic radio source. 

1830-210: 2EG J1834-2138 is identified with PKS 1830-210 with 98% confidence. No VLBI 
observation has been reported for this source. This 8 .Jy source has not been identified optically (it is 
6° from the Galactic disk near the Galactic center). A VLA radio image shows what appears to be 
a steep-spectrum Einstein Ring and two parity-reversed images of a flat-spectrum core-jet structure 
separated by 1" (Jaimccy et al. 1991). This structure is interpreted as a gravitationally lensed AGN 
jet. Thus, it could be an EGRET blazar. Wiklind and Combes (1996) established through mm 
observations of hydrocarbon absorption lines that the redshift of the lensing galaxy is z=0.89. If a 
detailed lens model can be developed, the autocorrelation function of future high-temporal density 
7-ray flux measurements may reveal information about the location of the 7-ray emitting region. The 
expected time delay between the two components is ^lO'' seconds (Jaunccy et al. 1991). 

1908-201: PKS 1908-201 is identified with 2EG J1911-1945 with 92% confidence. This 2 Jy source 
has not yet been identified optically (it is 21° from the Galactic center). The 8.4 GHz fiux density of 
PKS 1908—201 is 2.3 Jy (Wright et al. 1991). Using the consequent index a = (instead of assuming 
a = —0.5 without a 2.7 GHz detection) increases the probability of the identification to ~95%. The 
correlated VLBI 2.29 GHz flux density is 0.46±0.05 Jy (Preston et al. 1985). It would be worthwhile 
to observe the the parsec scale radio structure of this source. An optical (or IR) identiflcation and 
spectrum is also of interest. Preston et al. (1985) give a position which is accurate to 1": B1950, 19 8 
12.6 -20 11 57. or J2000, 19 11 09.8 -20 06 57. 

1933—400: This identification by TH95 is plausible, but narrowly fails to meet our criteria with 
p{id\r) 0.34 and Svlbi = 280 ± 30 mJy. It is possible that this 2.29 GHz VLBI fiux density 
measurement (Preston et al. 1985) is under representative because of variability. 

2022—077: This identification by TH95 is plausible, but fails to meet our criteria with p{id\r) = 0.61 
and Svlbi = 130 ± 20 mJy. It is possible that this 2.29 GHz VLBI fiux density measurement (Preston 
et al. 1985) is under representative because of variability. 

2155—304: An EGRET detection of this X-ray selected BL Lac type object was reported by 
Vestrand, Stacy and Sreekumar (1995). Since these authors do not provide a specification of the 95% 
confidence source location region, a likelihood analysis (Mattox et al. 1996a) of the VP 404 data for 
E> 100 MeV was done. The 95% confidence contour is well fit by an ellipse centered on 1=17.26 , 
b= —52.28, and with semi-major axis a=45', semi-minor axis b= 40', and 79°. The analysis of 
§2 indicates that B2155-304 (S'5=407 mJy, a = 0.3) can be identified with the EGRET source with 
p{id\r)=12% confidence. A test for variability of the EGRET source using the average flux of phase 
1&2 of the mission, VP 209, and VP 404 results in y^— 7.3 for 2 DOF. Thus, a constant 7-ray flux can 
only be rejected at the 97.5% confidence level. Preston et al. (1985) measured a 2.29 GHz VLBI flux 
density of 144±8 mJy. The identification is very plausible, but fails to meet our criteria for inclusion 
in Table 3. 

2200+420: A 4.4cr EGRET detection of 2200+420 (BL Lacertae) is reported by Catanesc et al. 
(1997). The 95% confidence contour for the event selection E> 100 MeV is well fit by an ellipse 
centered on 1=92.59 , b= —10.44, and with semi- major axis a=89', semi-minor axis b= 58', and (j)= 
35°. The analysis of §2 indicates that B2200+420 (^5=3.6 Jy, a = -0.2) can be identified with the 
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EGRET source with p{id\r)=99.3% confidence. A high confidence identification is obtained in spite 
of the large error region because the 5 GHz flux density is large. 



3.6 Limitations of this Analysis 

Although it is a significant improvement over previous EGRET identification work, we note that 
our technique is subject to a number of potential sources of systematic error. The error ellipses of TH95 
indicate statistical error only. Substantial systematic error is also expected, but was not given because of 
the difficulty of quantifying it. As discussed by TH95 and Mattox et al. (1996a), error in the model of the 

Galactic diffuse 7-ray emission and undetected nearby 7-ray point sources can have a substantial effect on 
EGRET position estimates. Because the present investigation is limited to |6| > 3°, the effect of error in the 
Galactic diffuse model of Hunter et al. (1996) is not severe for the EGRET sources we study in this paper. 

The radio flux densities of blazars vary by factors as large as ~10. Therefore, it is problematic to use the 
radio survey flux for the identification analysis because the flux at the survey epoch could differ significantly 
from the average flux. Also, the radio spectral index we use was obtained from non-simultaneous surveys at 
4.85 GHz and 1.4 GHz. It is subject to error because of flux variability. The error that this introduces can 
be mitigated by follow-up radio observations of sources which are potential EGRET source identifications. 
For this reason, potential identifications with values of p{id\r) as small as ~0.05 should not be dismissed 
until follow-up observations have been made. 

The necessary simplification of assuming that the a priori probability, ry, is independent of EGRET 
exposure and radio spectral index also contributes some error to our estimates of the probabilities of correct 
identifications, but we expect that they are not severely compromise. An additional minor error in our 
identification analysis may result from some deviation of the dependence of the log of the likelihood of the 
EGRET data on the source position from the assumed paraboloid. Also, a minor error in our identification 
analysis is introduced by the discernible but small differences between the elliptical fits and the actual 
position confidence contours of TH95 which are available on the ApJ CD ROM corresponding to TH95. 

For these reasons, the probabilities of correct identification derived here are indicative rather than 
definitive. We expect that for identifications with an a posteriori probability p(irf|r) > 0.9, l—p{id\r) could 
be in error by as much as a factor of 5. Thus, p{id\r) = 0.9 implies 0.5 < p{id\r) < 0.98. Similarly we 
expect that for identifications with an a posteriori probability p{id\r) < 0.1, p{id\r) could also be in error by 
as much as a factor of 5. Thus, p{id\r) = 0.1 implies 0.02 < p{id\r) < 0.5. Although the large uncertainty 
of p{id\r) is not ideal, it is not so large as to render this analysis useless. We expect that our method offers 
an order of magnitude improvement in precision over that of TH95. 

For latitudes less than ~15°, a substantial fraction of the radio sources will be Galactic radio sources. 
We search for counterparts in this range because correct blazar identifications could result nonetheless. In 
fact, one planetary nebulae and one H77 region appear in Table 1 (and are noted in the footnotes), but are 
not thought to be 7-ray source counterparts. We find that no excess of compelling identifications results 
at low latitudes. Therefore, undiscerned Galactic radio sources are apparently not a major problem for 
Galactic latitudes 3°< |&| <15°. 

We do not quantitatively incorporate in our identification analysis the variability of EGRET sources, 
nor the amount of correlated VLBI flux density — although these certainly pertain. They are considered in 
the selection of the high probability EGRET identifications of Table 3 and the "most interesting" potential 
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identifications of Table 4. 

Also, there is potentially more information for some sources which would be difficult to quantitatively 
incorporate in our identification analysis. We note that the analysis described in §2 results in a formal 
probability of 70% that MRK 421 is correctly identified as an EGRET source. This results from the 
low radio flux of MRK 421, and the large EGRET position uncertainty for this weak source. Additional 
information drastically improves the confidence of this identification. MRK 421 has distinguishing 
characteristics which effectively reduce the number density of such sources. In addition to a substantial 
correlated VLBI flux density and domination of the optical emission by continuum synchrotron radiation, 
MRK 421 shows a large and highly variable X-ray flux. It has been detected at TeV energies through 
the observation of Cerenkov light produced in the atmosphere by 7-ray induced electromagnetic cascades 
(Buckley et al. 1996). The TeV flux is highly variable and it has been observed to vary in correlation with 
the X-ray flux (Buckley et al. 1996). The TcV identification is very secure because the position of the TeV 
emission has been determined to ^4' and is consistent with the MRK 421 position (Buckley et al. 1996). 

The identifications of Table 1 indicate that the OVV property is common to many of the sources 
detected by EGRET which have been frequently observed with optical telescopes. Thus, the OVV property 
may serve to distinguish sources and improve identification confidence. This may also be true for optical 
polarization and BL Lac classification. However, careful work remains to be done to ascertain that these 
characteristics distinguish EGRET sources. Wagner (in preparation) finds that 80% of all sources which 
have a flat radio spectrum up to 90 GHz are OW irregardless of whether or not they are detected by 
EGRET. The fraction of EGRET sources which are believed to be OVV is not significantly larger than 
this; however not all have been yet been extensively monitored in the optical band. Also, the fraction of 
blazars detected by EGRET which are radio-selected BL Lac type objects is not significantly different from 
the fraction of flat-spectrum radio sources which are BL Lacs. Also, optical polarization may be generic for 
flat-spectrum radio sources. Fugmann (1988) has studied optical polarization observations of the Kiihr et al. 
(1981) sample. He finds that the probability of measuring significant optical polarization for a flat-spectrum 
radio source in a single measurement is 62%. It is possible that all flat-spectrum radio sources would have 
the property of high optical polarization in at least one observation if sufficient observations were made. 

4. Discussion 

We find that a substantial number of the identifications of TH95 are far from certain. This has 
implications for most studies of the "EGRET blazars" . It impacts the study of the 7-ray luminosity function 
and evolution (Chiang et al. 1995). It should be considered in studying the correlation of 7-ray and radio 
luminosity (Padovani et al. 1993; Stecker, Salamon, and Malkan 1993; Dondi and Ghisellini 1995; Miicke 
et al. 1996), and in the use of the EGRET detections in calculations of the production of isotropic diffuse 
high-energy 7-ray background by AGN (Padovani et al. 1993; Stecker et al. 1993; Salamon and Stecker 
1994; Chiang et al. 1995; Stecker and Salamon 1996; Kazanas and Pcrlman 1996). Careful consideration of 
the confidence of the identification is also in order when seeking distinguishing properties of the "EGRET 
blazars" . Unless explicit consideration is given to potential misidentification, we suggest that these analyses 
be restricted to the 42 EGRET blazars of Table 3 for which p{id\r) exceeds 90%, or a high confidence 
identification is otherwise indicated. We estimate that ^1 of these identifications is specious. Additional 
EGRET exposure and additional analysis and observation at other wavelengths may eventually contribute 
-^10 more sources to this population. 
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4.1 On the Detection of Blazars by EGRET 

We have demonstrated in §3.1 with confidence far beyond question that EGRET is detecting radio loud 
AGN. Six of the identifications in Table 1 have a probability of being correct of at least 99.9% (0420—014, 
0528+134, 3C 273, 3C 279, 1633+382, and 3C 454.3) All of these have flat radio spectra (a > -0.5). 
All have substantial correlated VLBI flux density. All 6 show significant changes in optical flux on short 
time-scales. Of these 6, the 4 which have been extensively observed show significant optical polarization. 
Also, the 5 which have been observed sufficiently with VLBI to discern it show superluminal motion of the 
parsec scale structure of radio jets. Therefore, it is apparent that EGRET is detecting the blazar class of 
AGN. 

The other identifications of Table 3 all have radio spectra with a > —0.5, except for 0521—365 for 
which a = —0.7 from the PMN catalog. This source has not been extensively observed. The report of an 
EGRET detection lead to VLBI observations which reveal that the spectrum of the core is 04.9-8.4 = +0.3 
(Tingay et al. 1996). Thus, if this source is truly stccp-spcctrum in single dish observations, it is due 
to the flux of the radio lobes. Therefore, it appears that all high confidence EGRET identifications are 
fiat-spectrum radio sources. Fugmann (1988) concludes from his study of the optical polarization that at 
least 2/3 of all flat-spectrum radio sources are HPQ sources and therefore blazars. Therefore, it is plausible 
that all extragalactic EGRET source are blazars. 

It is interesting that EGRET is detecting only ^10% of the blazars. It is possible that this results 
from time variability of the gamma-ray emission. However, a "duty cycle" of ~30% is observed for the 
strong EGRET sources — too large to explain the fact that only 10% of the blazars are detected, unless 
there is strong variability on time scales longer than ^1 year. Another possibility is that the gamma-ray 
luminous blazars differ intrinsically from the blazar parent population. E.g., the gamma-ray emission may 
be more narrowly beamed than the longer wavelength emissions so that the 7-ray blazar sub-population 
has a smaller angle between the jet and the line of sight than the blazar parent-population (Salamon & 
Stecker 1994, Dermer 1995). 

4.2 Non-Detection of Steep Radio Spectrum and Radio Quiet AGN by EGRET 

In Figure 6 we present a scatter plot of S5 and a for all radio sources considered. The identifications of 
Table 1 & 2 with p{id\r) > 0.02 are surrounded by circles. The size of the circle increases with the confidence 
of the identification. It appears that the distribution of observed spectral indices for the identified sources 
is consistent with that of flat-spectrum sources (a > —0.5) given the errors in the spectral indices due to 
the non-simultaneity of the radio surveys. Thus, it is plausible that EGRET is detecting only flat-spectrum 
radio sources. 

It is a possible concern that the non-detection of steep-spectrum sources results from our use of the 
reduced number density of flat-spectrum radio sources in the calculation of the likelihood ratio. Therefore, 
we have repeated the analysis described in §2 with a number density of radio sources piS^) which depends 
only on flux density. This results in values of LR smaller by a factor of ~2. However, significant detections 
persist and the list of strong identifications is nearly unchanged. The sum for the sources in Table 1 is then 
T,p{id\r) = 31.7 (down from 36.6 with the radio index dependent munbcr density described in section 2). 
In the mfll hypothesis, the Monte Carlo calciflation with random EGRET positions described in section 
3.1 yields Y,p{id\r) = 4± 2. Thus, the result that EGRET is detecting radio loud AGN does not depend 
upon the assumption that the flat spectrum sources are more likely to be detected. It follows that the 
steep-spectrum sources are generally not being detected. Therefore, the use of the spectral index in the 
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Fig. 5. — The distribution of the visibilities (VLBI flux density/total flux density). The 451 sources in the 
catalog of Preston et al. (1985) with a > —0.5 and a visibility are shown with the histogram using thin 
lines. Upper limits are treated as detections at the upper limit of the visibility. The 19 EGRET blazars with 
visibilities in the catalog of Preston et al. (1985) are shown with the histogram using bold lines. 
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Fig. 6. — The distribution of radio flux density {S5) and radio spectral index (a) for all sources in the PMN 
and GB surveys with 55 > 100 mJy. We plot an index of a = —0.5 if an index is not available from the 
survey catalogs (as assumed for the analysis of §2). The consequent enhanced source density at a = —0.5 is 
apparent. The identifications of Table 1 and Table 2 with p{id\r) > 0.02 are surrounded by circles (using a 
from Table 1 or 2). The size of the circle increases with p{id\r). The five different sizes correspond to the 
following ranges: 0.02 < p{id\r) < 0.1, 0.1 < p{id\r) < 0.5, 0.5 < p{id\r) < 0.9, 0.9 < p{id\r) < 0.99, 0.99 
< p{id\r) < 1.00. 
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calculation of confusing source number density is justified. 

We also suggest that EGRET is not detecting radio-quiet AGN. If an unidentified EGRET source 
is actually a blazar, that blazar is very likely to be among the possible identifications given here. If the 
EGRET flux density is variable, the likelihood of a blazar identification is enhanced. We note that every 
variable EGRET source in Table 1 (with variability index > 2 corresponding to confidence of variability 
>99%), has a plausible radio- loud identification. If we assume that the 7-ray emission of all AGNs would 
vary, we can conclude that EGRET is not detecting any radio-quiet AGN. We note that there are several 
unidentified EGRET sources at \b\ < 3° which are excluded from our analysis which are probably Galactic 
(e.g. Tavani et al. 1997). 

4.3 Suggested Observations at Other Wavelengths of Potential Gamma-Ray Selected Blazars 

Although we have found only one new EGRET identification in this analysis (PKS 1830—210) which 
qualifies for inclusion in the high confidence identifications of Table 3, it is apparent that a substantial 
number of the previously unidentified EGRET sources are also blazars. Further study at other wavelengths 
(primarily radio and optical) is required to ascertain whether or not some of these candidates have blazar 
properties. Most of these potential EGRET sources have not been studied extensively. Upon further 
investigation, they may be found to be "gamma-ray selected blazars". In Table 4 we list the 17 "most 
interesting" new potential identifications. This table also gives the best positions available for these sources. 

Further radio observation can establish variability and a radio spectrum from simultaneous observations. 
Radio interferometric images (VLBI, or VLA/AT) can be obtained to establish that the radio source is an 
AGN. A large correlated VLBI flux density is indicative of EGRET sources. Measurements of flux densities 

in the mm radio band may help to distinguish blazars. Measurements of optical variability and polarization 
can be used to establish membership in the OVV or HPQ blazar class. Spectra to establish redshifts, and 
equivalent widths of emission lines are of interest, as are X-ray and infrared observations. Simultaneous 
multiwavelength observation will be of interest once identifications are firmly established. In principle, the 
identification could be secured by observation of gamma-ray flux variation which correlates with variation 
at a longer wavelength. However, the EGRET sensitivity is generally not sufficient for this unless the 7-ray 
flux is exceptionally large. 

In many cases, an optical counterpart is obvious from an examination of the POSS (or ESO) plates for 
the GB or PMN error region. In other cases, an optical counterpart is not clearly apparent in this ^20" 
region because of source confusion. Radio observations to obtain an interferometric position with the VLA 
or the Austrahan Telescope Compact Array (AT) are then appropriate. Positions accurate to ~0.01" are 
given in Table 4 for sources of likely interest which were observed by Patnaik et al. (1992 & 1996) with the 
VLA for potential use as MERLIN phase calibrators. 

4.4 Pulsar Candidates 

The non-variable EGRET sources for which a radio-loud identiflcation is not apparent are more likely 
than other EGRET sources to be pulsars. Tables lfc2 may thus be used to to select optimal pulsar 
candidates. Romani and Yadigaroglu (1995) predict that 70% of gamma-ray luminous pulsars arc in fact 
radio-quiet. Direct detection of periodicity for more than ^2 through a FFT of the EGRET data is unlikely 
(Mattox et al. 1996b). More work is required to determine if the predicted spatial distribution of pulsars 
(Yadigaroglu and Romani 1996) matches that of the pulsar candidates of Tables 1&2. 

We note that both EGRET sources in Table 1 which coincide with SNRs (2EG J0008-F7307 and 2EG 
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J0618+2234) have no flat-spectrum radio identifications with p{id\r) > 2x 10~^. This is consistent with the 
claim (Stumer and Dermer 1994, Esposito et al. 1996, Sturner, Mattox, and Dermer 1996) of a correlation 

between EGRET sources and SNRs. This correlation could be either the result of enhanced interaction of 
cosmic rays with the interstellar medium near SNR, or pulsars associated with the SNRs or with young 
stellar associations co-located with a SNR. 

4.5 The Correlation of Gamma-ray Flux with Radio Flux 

With the limited EGRET data, the dependence of the a priori probability, rj, on ^5 cannot be 
understood in detail. A simple hypothesis is that ~20% of all flat spectrum radio sources emit 7-rays, and 
that the observed decrease in r] with radio flux is the result of detecting a fraction of 7-ray blazars which 
decreases with radio flux because 7-ray flux is correlated with radio flux and EGRET has a relatively high 
detection threshold. The initial reports of a correlation between 7-ray luminosity and radio luminosity 
(Padovani et al. 1993; Stecker et al. 1993; Salamon and Stecker 1994; Stecker and Salamon 1996) are flawed 
in that they fail to consider the effect of a common redshift on the apparent correlation. An analysis 
which overcomes this defect through partial correlation analysis (Dondi and Ghisellini 1995) indicated 
a correlation with 99.5% confidence. However, a partial correlation analysis which also considers upper 
limits (Miicke et al. 1996) does not find significant correlation. We present a simple argument below which 
establishes with 99.998% confidence that the maximum observed 7-ray flux of EGRET blazars correlates 
with the average 5 GHz radio flux. 

We simply note that all of the EGRET blazars which are 7-ray bright are also bright radio sources. 
The ten EGRET blazars with a peak 7-ray flux above lO^^^cm-^s-iat E>100 MeV (0208-512, CTA 26, 
0528-M34, 1156-^295, 3C 279, 1406-076, 1622-297, 1633-^382, 1730-130, and 3C 454.3) all have an 
average radio flux density S5 which exceeds 1.0 Jy. Prom the GB and PMN catalogs, we estimate that 58% 
of the 5*5 >1.0 Jy radio sources have flat spectra {a > —0.5), and that there are 183 of them away from 
the Galactic plane (|6| >3°). The fraction detected as bright 7-ray blazars is 0.05±0.02. Using Poisson 
statistics, a 99.998% confldence lower limit on the fraction of these which are bright 7-ray sources is 0.0098. 
Examining the GB and PMN catalogs over narrow intervals (250 < 55 < 350 mJy, 350 < S5 < 500 mJy, 
and 500 < 6*5 < 1000 mJy) and correcting for the variable fraction of sources which have measured spectral 
indices, and the variable fraction which are flat spectrum; we estimate that there arc 1130 flat spectrum 
radio sources with 250 mJy < ^5 < 1.0 Jy away from the Galactic plane (|&j >3°). If the 7-ray flux of 
EGRET sources were not correlated with radio flux, we would expect 5±2% or 62±20 of these to be bright 
7-ray sources. Using the 99.998% confidence lower limit on the fraction detected, at least 11.1 of these 
radio sources should be bright 7-ray sources. None are observed. The Poisson probability of observing 
when 11.1 are expected is 1.5xl0~^. We thus reject the hypothesis that the 7-ray flux of EGRET sources 
is not correlated with radio flux with 99.998% confidence. It is possible that the correlation is non-linear as 
suggested by Miicke et al. (1996). 

This analysis assumes that the distribution of a is the same for the sources where it was not measured 
as it is for sources where it was measured. This may not be true in detail, but it is expected to be a 
good approximation for sources with ^5 > 250 mJy. We have effectively analyzed the correlation between 
the maximum observed 7-ray flux of EGRET blazars and the average 5 GHz radio flux because that is 
essentially what the current data offers. Most EGRET blazars are only detected by EGRET during a 
flaring state. And the GB and PMN surveys with a 5 GHz observation at one epoch offer approximately an 
average flux. 

If bright EGRET blazars with 250 mJy < ^5 < 1.0 Jy existed, the method of §2 could identify them 
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because the 7-ray error regions would still be small. We analyzed the 10 bright EGRET sources listed 
above with the counterpart radio flux artificially altered to 55=250 mJy. All identifications remain above 

the threshold for inclusion in Table 1 (p{id\r) >0.GG1), and 9 of 10 above the threshold for variable EGRET 
sources in Table 4 {p{id\r) >0.02). Two sources (3C 279 and 0528+134) would still be above the threshold 
for variable EGRET sources in Table 3 {p{id\r) >0.3). This a mute point however. There are no bright 
EGRET sources (peak 7-ray flux above 10~'^cm-2s-\ E>100 MeV) at high Galactic latitude (|6| >3°) 
which do not coincide with either a pulsar or a radio source with ^5 >1.0 Jy. In fact, this argument can 
be extended to the myriad radio sources with 5*5 < 250 mJy. The prevalent interpretation of the isotropic 
high-energy 7-ray diffuse radiation is that it is due to the 7-ray emission of these sources (Padovani et 
al. 1993; Steckcr et al. 1993; Salamon and Stecker 1994; Chiang et al. 1995; Stecker and Salamon 1996; 
Kazanas and Perlman 1996). If this interpretation is correct, there must be a correlation between 7-ray and 
radio flux or else the weakest radio sources would also be detected as bright 7-ray sources. 

We note that a Spearman rank-order correlation test of the sources in Table 3 for the 5 GHz radio 
flux of Table 3 and the peak 7-ray flux (without doing a K-correction) suggests a correlation (with 67% 
confidence). The same test shows a correlation with 99.2% confidence for the peak 7-ray flux and correlated 
VLBI flux of Table 3. This is consistent with the finding of §3.3.1 that the EGRET blazars are distinguished 
by a large VLBI flux density; and adds more weight to arguments that the 7-ray emission is taking place at 
the base of a relativistic jet. A tighter correlation might be obtained using simultaneous 7-ray and VLBI 
fluxes (perhaps allowing for a delay in VLBI flux because of synchrotron self-absorption of radio flux at the 
base of the jet). We note that the same analysis with the non-simultaneous 8.4 GHz VLBI flux densities of 
the 28 sources in Table 3 which were measured by Morabito et al. (1986) shows much less correlation (the 
Spearman rank-order correlation test indicates a 34% probability of a correlation). To some extent this is 
caused by the reduced number of sources included in the correlation analysis. Perhaps it also reflects higher 
radio variability at higher frequency. 

5. Summary 

We present a means to quantitatively assess the probability that an identiflcation of an EGRET source 
with a radio source is correct. We also demonstrate conclusively that EGRET is detecting 7-ray emission 
from the blazar class of AGN. Our analysis suggests that EGRET is not detecting radio quiet and steep 
radio spectrum AGN. This identification method is used to ascertain with 99.998% confidence that the 
7-ray and radio fiux of EGRET blazars are correlated. Previous EGRET identifications are re-examined. 
Table 3 gives the 42 EGRET blazars for which the confidence of a correct identification is high. It is found 
that the peak 7-ray fiux and the VLBI radio flux of these sources is correlated with 99% confidence. One 
high confidence AGN identification (PKS 1830—210) is found among the unidentified EGRET sources of 
TH95. Sixteen additional potential identifications are presented in Table 4 with the best available position 
for study at other wavelengths. 
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